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ABSTRACT 

Aims. In order to directly study the role of magnetic fields in the immediate vicinity of protostars, we use Very-Long-Baseline 
Interferometry (VLBI), aiming at the detection of non-thermal centimetric radio emission. This is technically the only possibility to 
study coronal emission at sub-AU resolution. 

Methods. We performed VLBI observations of the four nearby protostars HL Tau, LDN 1551 IRS5, EC 95, and YLW 15 in order to 
look for compact non-thermal centimetric radio emission. For maximum sensitivity, we used the High Sensitivity Array (HSA) where 
possible, involving the Very Long Baseline Array (VLBA), the phased Very Large Array (VLA), as well as the Arecibo, Green Bank, 
and Effelsberg radio telescopes. 

Results. While all four protostars were detected in VLA-only data, only one source (YLW 15 VLA 2) was detected in the VLBI 
data. The possibility of non-detections due to free-free absorption, possibly depending on source geometry, is considered. For 
YLW 15 VLA 2, the prospects for an accurate orbit determination appear to be good. 

Key words. Stars: pre-main sequence, Stars: coronae, Stars: magnetic fields, Radio continuum: Stars, Stars: individual: YLW 15 



1. Introduction 



Accor ding to the classification scheme suggested by lLadal 
(fl987h . low -mass pre-main sequence (PMS) stars form an evo- 
lutionary sequence. The seq uence starts with t he youngest (age 
< 10 4 y) "class 0" sources dAndre et alJI 1993b . that are not de- 
tected at near-infrared wavelengths and are deeply embedded 
in very dense dust envelopes of low temperature. It then pro- 
ceeds over less embedded class I sources (age 10 5 y), to virtu- 
ally envelope-free, optically visible, class II protostars (classical 
T Tau stars, 10 6 - 10 7 y). After being class III protostars or weak- 
lined or naked T Tau stars (i.e., without or with only very little 
disk material), they reach the zero-age main sequence at ages 
> 10 7 y. For a discussio n of radio pro perties of Young Stellar 
Objects (YSOs), see, e.g jAndrel d 1996b . 

Powerful X-ray and gyrosynchrotron radio-continuum flares, 
orders of magnitudes stronger than solar flares, indicate that 
protostars exhibit unusually h igh levels of magnetic activity. 
Feig elson & Montmerld ( 1 1999b identify four possible magnetic 
geometries in which hard X-ray and non-thermal radio emission 
might arise: 

(1) Solar-type multipolar fields with both foot-prints rooted in 
the stellar photosphere. 

(2) Field lines connecting the star to the circumstellar disk at 
its co-rotation radius. This could be the explanation for the 
quasi-periodic X-ray flares observed towards the class I pr o- 
tostar YLW 15 dTsuboi et al.ll2000t iMontmerle et al.ll2000b . 

(3) Field lines above the co-rotation radius. 

(4) Magnetic loops with both origins in the disk. 
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The presence of a disk is fundamental for magnetic energy 
release in the three last cases, but irrelevant for the first 
case. Further possible configurations include interacting coronal 
structures in binaries (see below). Apart from gyrosynchrotron 
emission, also completely different non-thermal emission mech- 
anisms can occu r, e.g.. coherent emission like electron cyclotron 
maser radiation dDuIkJl 1985t ISmith et al.H2003l) . 

In the late 1980s, the first VLBI observations of YSOs 
were reported. iFelh et all dl989h described the VLBI detection 
of radio emission from the 6 l Ori A multiple system; 8 l Ori A 
is one of the "Trape zium stars" at the center of the young 
Orion Nebula Cluster. Andr e et al ] d!99ll) reported non-thermal 
gyrosynch rotron radiat i on tow ards the young magnetic B star 
p Oph Sl. lAndreetaTJ dl992l) published a VLBI survey of the 
p Oph star-forming region, showing that several T Tauri stars 
have non-thermal radio emission, although detections did not in- 
clude embedded protostellar sources (like IRAS 16293-2422). 
The very active system HD 283447 (= V773 Tau), located in 
the Taurus-Auriga molecular cloud compl ex at a distance of 
d- 14 pc, could be resolved at VLBI scales (Ph illips et al.ll99ll 
1996). It consists of two class III objects. The VLBI observations 
have shown that the emission arises from regions of some tenths 
of an AU, and that the bulk of the radiation comes from two 
small unresolved components related to the two stars. The rea- 
son of the large flares observed in this system is different from 
the mechanisms listed above; the flares are clustered around the 
periastron passage of the binary system and therefore probably 
due to colliding in terbinary stellar corona l structure involving 
helmet stream ers (Mass retaill2002l [20 06). As a matter of fact, 
iPhilhps et all dl99llll996h ~"have performed several other VLBI 
observations of this source outside the periastron flaring window 
and observed a source of a size of a few stellar radii, or anyway 
much smaller than the dimensions of the binary system. 
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In order to study the predicted role of the disk as described 
by iFeigelson & Montmerld d!999l) . it is perhaps most interest- 
ing to observe younger objects than class III. This is difficult, 
since only very few candidates are known (i.e., protostars show- 
ing both X-ray and non-thermal radio emission). It is also possi- 
ble that a putative weak radio source is embedded in an absorb- 
ing plasma showing thermal (free-free) emission. The magnetic 
fields which lead to magnetic reconnection and generate flares at 
small scales (< 0. 1 AU) play an important role in driving thermal 
radio jets and bipolar outflows; the shock front produced by the 
collision of the outflows with the ISM ionizes part of the gas and 
creat es ionized regions app earing elongated at VLA resolutions 
(e.g.. lTsuiimoto et aT]|2004l) . Any emission process taking place 
at smaller scales will suffer from significant free-free absorption 
and will be difficult to detec t. This is exactly the mechanism in- 
voked by Gira rt et ail d2004l) to explain their non- detection with 
the VLBA at 6 cm of the class I protostar YLW 15. lAndreldl987l) 
estimates how easily emission from within this region can be en- 
tirely concealed by ionized material: A mass-loss wind of only 
10" 11 M yr~' suffices to entirely absorb 5 mJy of flux density 
of T/, > 10 7 K from behind. In conclusion, the real problem in 
VLBI observations of protostars younger than class III is that the 
weak non-thermal emission suffers from strong free-free absorp- 
tion. Additionally, free-free emission observed at angular scales 
probed by the VLA may be "resolved out" at VLBI scales. The 
shortest baseline in the VLBI data discussed here corresponds to 
angular scales of about 07 1 . Thus, already emission of slightly 
larger extent remains undetected. However, the rare detections of 
non-thermal radio emission towards protostars may also reflect 
an intrinsic property of these sources. 

The VLBI observations (using the VLBA, the VLA, and 
the Effelsbe r g 100 -m telescope) at a wavelength of 3.6 cm by 
Smit h et all (|2003) of the class II protostar T Tau-S are very 
important in this respect: Weak and variable emission (between 
1.4 and 3.8 mJy) was finally detected, even if unresolved with 
the 1 mas HPW beam. From multi-epo ch VLBI observations of 
the same source. lLoinard et"aT] ((2005) confirm the detection of 
a compact radio source and derive a parallax yielding a precise 
distance (141.5^ pc). 

In this paper, we report the results of our search for compact, 
coronal radio emission towards protostars, using high-sensitivity 
VLBI observations. In Sect. [2] we describe our target selection 
procedure before a description of the actual VLBI observations 
and their analysis in Sect. [3] The results are described in Sect. [4] 
and finally, we summarize our findings in Sect. [5] 



2. Target selection 

In order to establish the role of the disk or star/disk interaction, 
we planned to observe nearby protostars at an earlier stage than 
class III at A — 3.6 cm using VLBI, trying to detect non-thermal 
radio emission. Given the difficulties described, we first had to 
define a good sample of target sources. We selected class I and 
class II protostellar sources which were detected in centimet- 
ric radio as well as in (variable) X-ray emission as tracers of 
coronal magnetic activity. Since information on polarized radio 
emission from protostars — a sign for non-thermal emission — is 
sparse, sources were selected that had either variable radio emis- 
sion or a negative spectral index indicative of non-thermal radio 
emission. Only few sources are known to fulfil these criteria so 
that no additional criteria (e.g., distance) were applied. In our ac- 
tual observations (see below), synthesized beam sizes of down to 
0.2 x 0. 1 AU were realized. We chose an observing frequency of 
8.4 GHz (A = 3.6 cm) as a compromise between a reduced influ- 



ence of free-free-absorption (compared to longer wavelengths) 
and ensuring a good sensitivity by using a standard frequency 
allowing the use of the Arecibo 305m telescope. 

Our sample consists of the class 0/1 binary protostar 
YLW 15, the class I binary protostar LDN 1551 IRS5, the class II 
source HL Tau, as well as the proto-Herbig Ae/Be star EC 95. 
Distances range from 130 pc to 310 pc (see Table [T]) . As noted 
above , YLW 15 was observed with VLBI before dGirart et all 
2004). Except for this source (YLW 15), all sources are in the 
declination range accessible to the Arecibo 305m telescope and 
thus were observed with the full High Sensitivity Array (HSA), 
consisting of the VLBA, the phased VLA, Arecibo, GBT, and 
the Effelsberg 100-m telescope. The radio luminosities of these 
sources, as deduced from previous VLA observations, are all 
very similar, as will become clear in the source descriptions. 



2. 1 . The binary protostar YLW 15 

YLW 15 is a binary protostar whose components are in the 
class 0/1 stages. It is located in the p Oph dark cloud at a dis- 
tance of 130 pc dRebull et al.ll2.Q04l) . and it shows X-ray and ra- 
dio emission. It was first found to be a radio binary ( with com- 
ponents named VLA 1 and VLA 2, iGirart et ai1l2000l), whose 
orbital parameters were estimated by ICuriel et al.l ([2003 ) from 
multi-epoch VLA observations. After an early VLBI attempt 
to detect non-thermal radio emission from t he YLW 15 system 
yielded an upper flux density limit of 1 mJy ( Andre et alii 19921) . 
recent VLBI observations were reported by IGirart et al.l (2004) 
who observed at a wavelength of 6 cm and achieved an upper 
limit of 0.2 mJy (4<x). However, the authors note that both the 
non-detection of the source on the VLA-Goldstone baseline and 
their finding that the source has a non-variable radio spectrum 
point towards an extended source of thermal free-free emission. 
Interestingly, their additional near- and mid-infrared data indi- 
cates that VLA 1 is more deeply emb edded (or less lum inous) 
than VLA 2 while the orbital motion dCuriel et alj|2003l) indi- 
cates that VLA 1 is the more massive object. 

VLA 2 is responsible for strong quasi-periodic X-ray emis- 
sion explained as due to the fast rotation of the star with re - 
spect to the disk dTsuboi et al.l 120001: iMontmerle et all 12000). 
The consequent star-disk shearing of the magnetic field lines 
gives rise to magnetic reconnection and therefore creation of 
energetic electrons producing the quasi-periodic X-ray emis- 
sion. Non-thermal radio emission is therefore expected. The in- 
tegrated flux at 3.6 cm, d erived from VLA measurements, is 
S3.6cm = 0.64 . . .0.78 mJy dGirart et al.ll2004 . When comparing 
the VLBA-only non-detection at a wavelength of 6 cm with VLA 
data at the same wavelength, it seems that at least two thirds of 
the total emission is thermal emission from a collimat ed outflow 
which is indeed partially resolved at VLA scales dGirart et al.l 
2004). The non-thermal contribution is thus embedded in an ab- 
sorbing thermal medium. 



2.2. The binary protostar LDN 1551 IRS 5 

This source, located at a distance of 140 pc in Taurus, is one of 
the best known deeply embedded young stellar objects. It drives 
one of the most spectacu lar bipolar outflows, the first ever de- 
tected dSnellet alJll9 80). extend ing over more than 10' (e.g. 



iMoriartv-Schieven & Snelll988h . The dark cloud LDN 1551 it- 
self is one of the nearest and most active regions of ongoing low- 
mass stars formation. IRS 5 is the most luminous YSO in this 
region (L = 40 L ), and deeply embedded (Ay ~ 150 mag, e.g., 
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Table 1. Journal of observations 1 



Source 


d[pc] 


Phase ref . 


Project 


Date 


Time (UTC) 


On-source 2 


VLA J 


Comments 


YLW 15 


130 


J1625-2517 


BF083 


30 Apr 2005 


06:15-10:15 


142.0 min 


B 


VLBA without LA 4 


HL Tau & LDN 1551 s 


140 


J043 1 + 1731 


BF084a 


03 Jul 2005 


12:30-16:30 


64.1 min 


BnC 


Arecibo data lost 


EC 95 


310 


J1838+0404 


BF084b 


29 Jul 2005 


01:00-04:00 


90.4 min 


C 




HL Tau & LDN 155 I s 

1 :„j ...:.u .u„ tt 


140 


J043 1 + 1731 


BF089 


24 May 2006 


15:00-19:00 


65.5 min 


AnB 


BF84a reobserved 



1 carried out with the High Sensitivity Array (HSA) except for BF83 (VLBA+VLA+GBT) 

2 scan minutes, not taking into account the number of antennas observing 

3 configuration of the VLA during the observation 

4 the Los Alamos VLBA antenna 

5 IRS 5 



IWhite et al.l2000l) . lRodrfguez et al.ldl998l) disc overed that IRS 5 
is a p rotobinary separated by 40 AU (see also [Rodrigue z et alj 
1986). Based on VLA observations at A — 7 mm, they detected 
two protoplanetary disks with semimajor axes of about 10 AU. 
The integrated flux densities at A = 3.6 cm are 0.78 mJy and 
0.70 mJy for the northe rn and southern component, respectively 
dRodrfguez et alj|l998l) . Most recently, Eim & Takakuwal (2006) 
analyze the source at 7 mm using the VLA combined with the 
Pie Town VLBA antenna, and find that LDN 1551 IRS 5 is actu- 
ally a triple system. 

iFavata et al.l d2002l) detected a faint X-ray source near IRS 5 
in an observation carried out with XMM-Newton. Ball y et aTl 
(2003), presenting Chandra X-ray data of higher angular resolu- 
tion (but still not resolving the binary), conclude from a position 
offset and a difference in absorbing column density that an X- 
ray source is located at the base of the HH 154 jet. The source, 
named ACIS 31, is located about Of.' 5 west of IRS 5 with 60 
counts registered in the 0.5 - 8 keV band in a 79 ksec expo- 
sure, 47 of them coming from within 1'.'7 of the source. Due to 
the difference in angular resolution between the XMM-Newton 
and Chandra data, it cannot be excluded that the X-ray emission 
is instead due t o the protostar itse lf (although that appears to 
be unlikely; see lFavata et al.ll2006l who interpret morphological 
changes in the source as due to jet motion). 

2.3. The class l/ll transition protostar HL Tau 



HL Tau is a flat (optical — far-infrared) spectrum T Tauri star, 
deeply embedded in a dusty nebulosity associated with a re- 
flection nebula. In the evolutionary sequence, it is in a transi- 
tion phase between class I and class II. There is strong extinc- 
tion t owards the central source (A v > 22 mag, Stapelfeld t et al.1 
1 1995b. although the star is detected in th e near-infrared 
dWeintraubetal.il 19951; iBeckwith & Birklll995l) . The morphol- 
ogy of the optical and NIR reflection nebula as well as the 
flat spectral energy distribution of HL Tau can be explained by 
scattering and thermal emission from an infalling, dusty proto- 
stellar envelope typical of Taurus class I sources ( Calv etet al.l 
Il994l) . The star is associated with the outflow system HH 266. 
At centimeter radio waveleng ths, HL Tau has been analyzed by 
iRodrfguez etakl (1 19921 11994 . using the VLA. At A - 1.3 cm, 
they observe a structure elongated approximately in the north- 
south direction surrounding the star (about 1.4 x 0.6 arcsec, or 
about 200 x 80 AU at a distance of 140 pc). The authors ar- 
gue that the size and the integrated flux density of this source 
(2.9 ± 0.4 mJy) at A = 1.3 cm indicate that the emission orig- 
inates from dust in a protoplanetary disk. At A = 3.6 cm, the 
picture is strikingly different: Here, the observed structure is 
elongated approximately in the east-west direction (with an in- 
tegrated flux density of 0.52 + 0.02 mJy). The authors argue 



that in contrast to the emission at A = 1.3 cm, the emission at 
A = 3.6 cm is dominated by free-free radiation from an ion- 
ized outflow. In order to study the structure of its circumstellar 
material, HL Tau has been s tudied extensively at millimeter and 
submilhmeter wavelengths dLav et al. 1994; Mu ndv et"ai] [l996; 
lLav et al.ll 19971; lLoonev et al.ll2000l) . lFavata et al.l d2003l) discuss 
X-ray data of HL Tau. While the source seems to have been 
quiescent during XMM-Newton observations, there are Chandra 
data containing evidence for a short-duration flare. 



2.4. The Proto-Herbig Ae/Be star EC 95 

EC 95 is a very young intermediate-mass star (i.e., up to 5 M ), 
presumably a precursor of a Herbig AeBe (HAeBe) star, there- 
fore not a low-mass protostar like our other sources. It is located 
in the Serpens cloud core at a distance of 310 pc. 

The source has an X-ray luminosity of 3.7 x 10 31 ergs -1 



(Preibisch 2003). Radio emission has been detected (with peak 
flux densities S%a ghz ~ 0.6 mJy and Ssghz ~ 0.8 mJy) by 
Smit h et ail {1999). The slightly negative spectral index (in S oc 
v") is inconsistent with thermal emission fr om an optically thick 
wind, where rather a — +2 is expected (Simon et "all 11983b . 
There is also evi dence for variability on timescales of years. 
Eiro a et all d2005l) report a peak flux density of 1.65 + 0.1 mJy 
on 24 Oct 1993, also from 3.5 cm observations carried out with 
the VLA in CnD configuratio n. Only about one month earlier (in 
the observations reported by Smi th et alii 19991) . the source was 
more than three times fainter. 

Although intermediate-mass stars are thought to be fully ra- 
diative, and therefore should not have a corona, previous de- 
tections of X-rays from HAeBe stars have been m ade (e.g., 
IZinnecker & Preibisch|[l99ll) . iTout & Pringlel (1 19951) suggested 
that radially differential rotation in intermediate-mass stars could 
lead to the formation of a temporary convective layer and 
hence a corona. They estimated that the lifetime of such a 
corona would be of order 10 6 years, approximately the timescale 
for intermediate-mass objects to evolve onto the main se- 
quence. On the other hand, a companion star responsible for 
the corona cannot be excluded. Another explanation of X-ray 
emission from Herbig A e/Be stars involves star-disk interactions 
(lHamaguchi et al.ll2005l) . 



3. Observations and data analysis 

In order to reach the highest possible sensitivity we used the 
world's biggest radio telescopes together with the NRAO Very 
Long Baseline Array (VLBA, ten antennas of 25 m diameter 
each), most notably the so-called High Sensitivity Array (HSA) 
consisting of the Arecibo 305 m, the Effelsberg and Green Bank 
(GBT) 100-m radio telescopes, the phased Very Large Array 
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Fig. 1. VLA-only maps of the four program sources YLW 15-VLA2, EC95, HL Tau, and LDN 1551 IRS5. Note that YLW 15 is 
dominated by VLA 1 (see also Fig. [2] for an image restored with a circular beam smaller by about a factor of two). For the two 
components of YLW 15, the proper motions as read from Fig. 2 in lCuriel et alJ (120031) are shown ('x' for VLA 1 and '+' for VLA 2; 
years 1990, 2000, and 2002 from north to south) together with the correlation position ('□'), see text. For all sources the contour 
lines delineate multiples of the 3cr image noise, (-2,- V2,-l,l, V2,2,...) x 3<x, increasing by factors of y2. The respective 3<x 
levels are 72 /jjy (YLW 15), 79 fjy (EC 95), 72 /jjy (HL Tau), and 74 fjy (LDN 1551). 



(VLA, equivalent in collecting area to a 130 m dish) as well 
as the VLBA. The usage of the phased VLA, equivalent in sen- 
sitivity to a large single-dish telescope of 130 m diameter, has 
the additional advantage of providing regular VLA data as well. 
Thus, the source positions can be checked before correlation. 
All runs were phase-referencing observations in dual polariza- 
tion (RCP & LCP) taken at 256 Mbit s" 1 , eight baseband chan- 
nels of 8 MHz bandwidth covering an aggregate bandwidth of 
32 MHz, and 2-bit sampling (mode 4cm 256-8-2-UL). Each ob- 
servation was correlated with 16 channels per IF. The field of 
view is limited by time-smearing due to the correlator averag- 



ing time. Its size is estimated to be ~ l'/8 for YLW 15 and 
~ 3'.'3 otherwise. Basic observation parameters including the 
phase reference sources used are listed in Table Q] The phase 
reference sources were at 52' (YLW 15), ~ 40' (HL Tau & 
LDN 1551 IRS5), and 3.6° (EC 95) of the program sources. We 
spent about three minutes on the program source and 1 .5 minutes 
on the reference source in every cycle. 

In order to obtain accurate positions for the correlation of 
the VLBI data, the VLA data for all sources were analyzed first 
(Fig. |2j see text below). Confirming the position was especially 
important in the case of YLW 15 due to the apparent orbital 
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Source 


RA (J2000) 


Dec (J2000) 


coravg 1 


Source 


synth. beam 




map peak 1 


map rms 1 


YLW15 VLA2 


16h27m26.89s 


-24°40'50.9" 


4s 


YLW 15 


1'.'3 x a/6, PA 


-10°3 


1.48 mJy- 


0.024 mJy 


EC 95 


18h29m57.92s 


+01°12'46.0" 


2s 


EC 95 


3'/ 1 x 2/4, PA 


-38°5 


0.55 mJy 


0.026 mJy 


HLTau 


04h31m38.432s 


+ 18°13'57.485" 


2s 


HLTau 


a/9 x a/3, PA 


-74°2 


0.40 mJy 


0.024 mJy 


LDN1551 2 

1 


04h31m34.158s 


+ 18°08'04.760" 


2s 


LDN1551 3 

— I : — tz r 


a/9 x a/3, PA 


-74° 2 


0.99 mJy 4 


0.025 mJy 



IRS 5 



motion of its components. All VLA and VLBI data were ana- 
lyzed with the NRAO Astronomical Image Processing System 
(AIPS). The VLA data of YLW 15 and EC 95 were amplitude- 
calibrated using a standard total flux density of 3C 286 of 5.22 Jy 
while HL Tau and LDN 1551 IRS 5 were calibrated using a stan- 
dard total flux density of 3C 84 of 3.17 Jy. While EC 95 and 
HL Tau were correlated simply at the positions of the singular 
emission peaks, and LDN 1551 IRS 5 was correlated at the posi- 
tion of the northern component, YLW 15-VLA2 was correlated 
at a position in the south-western extension to the main source as 
seen more clearly in a VLA-only image restored with a circular 
(y.'4 beam (Fig.|2]i. The correlation parameters used are listed in 
Table[2]for all sources. 

The VLBI data were analyzed as follows. After correcting 
the antenna parallactic angles recorded in horizontal coordinates 
and after applying the latest Earth Orientation Parameters, the 
dispersive delays caused by electrons in the earth's ionosphere 
were accounted for using Global Positioning System (GPS) 
models of the electron content of the ionosphere. For amplitude 
calibration, the T sys and gain curves of the phased VLA (which 
is recording ratios of antenna to system temperatures) were read, 
and the flux density scale was set using a calibrator source. Data 
taken during telescope slewing were flagged. After correcting 
sampler voltage offsets using autocorrelation total-power spec- 
tra, the a priori amplitude calibration was applied based on the 
r sys and gain tables. In order to remove instrumental residual de- 
lays, phase corrections had to be applied, ideally using pulse-cal 
information inserted into the data during the observations. Since 
there is no such information in the VLA and Arecibo data, these 
corrections had to be derived from the fringes on a strong source 
for baselines involving these telescopes. From them, the resid- 
ual delay and phase can be computed which can then be taken 
into account. Finally, global time-dependent phase errors have to 
be corrected (e.g., phase gradients between the IFs). These are 
mainly due to inaccuracies in the geometrical delays assumed 
during correlation, causing phases varying rapidly with time. 
They were corrected with the fringe-fitting algorithm devised by 
ISchwab & Cotton! (11983). Before the final application of calibra- 
tion information, the program sources were phase-referenced to 
their respective calibrator sources. As a check of global calibra- 
tion quality, the calibrator sources were referenced to themselves 
and subsequently imaged. No polarization calibration was at- 
tempted because of the very weak flux densities detected. Maps 
were created without a uv taper, setting the ROBUST parameter 
in IMAGR to zero. 



4. Results 

4.1. VLA data 

As stated above, we analyzed the VLA-only data taken of each 
source as part of the VLBI experiments in order to refine the 
source positions for the correlator. All sources were detected 
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16 27 26.96 26.94 26.92 26.90 26.88 26.86 

RIGHT ASCENSION (J2000) 

Fig. 2. VLA-only map of YLW 15, restored with a circular CC.'4 
beam. For the two com ponents of YLW 15 , the proper motions 
as read from Fig. 2 in ICuriel et al.l d2003l) are shown ('x' for 
VLA1 and '+' for VLA 2; years 1990, 2000, and 2002 from 
north to south) together with the correlation position ('□'). The 
contour lines delineate multiples of 0. 1 mJy ( ~ Act), increasing 
by factors of V2 (as explained for Fig. [T). 



(Table [3]). The VLA results for the program sources are shown 
in Fig.[TJand summarized in Table[3] 



YLW 15 For YLW 15 with its previously analyzed proper mo- 
tion, it was especially difficult to extract the position of YLW15- 
VLA2 out of the VLA data where VLA1 and VLA2 are blended. 
We chose the south-western extension to the main source VLA1, 
as seen more clearly in an image of our VLA data restored with a 
circular 0'.'4 beam, smaller than the synthesized beam (Fig.|2|i. It 
remains unclear whether this emission is du e to VLA 2. This p o- 
sition is unexpected given the predictions of lCuriel et al.ll2003l — 
the positions from their Fig. 2 are plotted into our VLA maps 
in Fig. [TJ and Fig. [2] — but the angular resolution of our VLA 
data is not good enough to decide this issue. However, already 
from the previous data it becomes clear that VLA 2 was mov- 
ing westward relative to VLA 1, and it is west of these previ- 
ous positions where we find emission. On the other hand, the 
predictions are based on only three positions, thus the orbit is 
no t yet fully de t ermin ed. Given that the observat ions discussed 
by ICuriel etail d2003l see also iGirart et al.ll2004l) were carried 
out in 1990, 2000, and 2002, the big change seen since then 
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in only about three years — if VLA2 is in the above-mentioned 
extension — could possibly be explained by the periastron pas- 
sage of VLA 2. More observations are needed in order to reliably 
determine the orbit. 

The integrated flux densities of th e two components (VLA 1 
and 2) as given bv lGirart et al.l d2004l) are 2.18 mJy and 2.14 mJy 
in 2000 and 2002, respectively; whereas the integrated flux den- 
sity of VLA 1 is 1.40 mJy and 1.51 mJy for the two epochs. In 
our case, the integrated flux density for both sources is 2. 1 8 mJy. 
Assuming therefore no significant long-term variability, the flux 
density of VLA 2 is 0.67. . . 0.78 mJy. 

LDN 1551 IRS 5 Both components of LDN 1551 IRS 5 are 
of similar flux density. We found the northern component, of 
special interest here, to be slight ly brighter than the so uthern 
component in our VLA data while Rodriguez et al. ( 1998) found 
the opposite in 1996 data. The peak flux density in our data of 
0.99 mJy at the position of the northern component compares to 
a peak flux density of ~0.5 mJy at the position of the sou thern 
component in their data (Fig. 2 of Rodrig uez et al.| [l998i). The 
source morphology in both datasets is similar, and differences in 
the observed relative source positions are much smaller than the 
beam size. For these reasons and because of the very elongated 
beam shap e, we refrain from disc uss ing proper motions in th is 
system (cf . iRodrfguez et al.ll2003l and lLim & Takakuwall2006l) . 

HL Tau The peak flux density of HL Tau, as determined from 
our VLA-only data is 0.40 mJy which is com patible with the 
peak flux density of ~ 0.3 mJy determined by Rodriguez et al.l 
( 1994, Fig. 2). Since their measurement in 1992, the source has 
shifted eastward by about 075. However, this displacement is 
along the major axis of the elongated synthesized beam in our 
data which also prevents a confirmation of the slight east-west 
elongation of HL Tau apparent in the 1992 data. 



EC 95 In 1995, EC 95 was observed twice at peak flux densities 
of 0.64 + 0.02 mJy (01 May, D array) and 0.62 ± 0.02 mJy (20 
Jul, A array), respectively; an observation carried out on 21 Sep 
1993 (in Cn D array configurat ion) yields a peak flux density of 
0.5 ±0.1 mJy dSmith et al.ll999l) . These flux density is very close 
to the 0.55 ± 0.07 mJy found in our data, alth ough different from 
the 1.65 ± 0.1 mJy reported by lEiroa etaTI (120051) for October 
1993. 



4.2. VLBIdata 

While all four target sources were detected with the VLA (ob- 
serving in various configurations, see above), the clean maps 
produced from the VLBI data contain only one detection. Basic 
results for each source, including the rms noise level achieved, 
are listed in Tab. |4] 

In the VLBI data of YLW 15-VLA2, we find a persistently 
detected weak source (see Fig. [3). In particular, the detection 
is reproduced when using different cell sizes in image recon- 
struction and different weighting parameters in the uv plane 
(ROBUST parameter in AIPS or using a uv taper to reduce the 
weights of the visibilities with the longest baselines), yielding a 
significance of 7 . . . 9<x. 

The peak flux density of this unresolved source is 0.145 ± 
0.016 mJy, therefore almost 20% of the flux density of VLA 2 is 
due to compact emission of < 0.4 x 0.1 AU (FWHM beam size). 
The flux density is too low for a discussion of detections on sin- 
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Fig. 3. Weak source found in the VLBI data close to the puta- 
tive position of YLW 15 YLW 2 at RA 16h27m26.88875s DEC 
-24°40'50.937". The source reaches a significance of ~ 9cr. The 
contour lines delineate multiples of 3<x, increasing by factors of 
V2 as explained for Fig. Q] The lcrrms noise levelis 15.4/Jy. 



gle baselines, but the use of all phase-referenced data on all base- 
lines for the full time of the observation allows the detection of 
the source. Towards VLA 1, no compact emission was detected. 
Interestingly, the direction of the outflow towards YLW 15 as 
seen in CO appears to be quite close to the direction of the line 
of sight be cause the integrated lin e wing emission overlaps con- 
siderably (Bontem ps et al.|[l996l) . This may cause a favourable 
viewing geometry onto the central source along the jet. 

YLW 15 has been observed in the near-infrared using ISAAC 
at the ESO Very Large Telescope (VLT) only 1.5 months after 
our VLBI run (16 June 2 005). Since only VLA 2 is supposedl y 
detected in J band (e.g., lAllen et all [20021: iGirart et alj|2004l) . 
these data can help assess the position of this source indepen- 
dently. Six exposures of 15 sec duration were obtained from the 
ESO Science Archive (program code 075.C-0561, PI: B. Nisini). 
They were reduced using the jitter algorithm implemented in the 
ESO Eclipse package. The accuracy of the astrometry was im- 
proved using a fit to the five 2MASS PSC sources in the field 
of view. With such a low number of reference sources, however, 
it is difficult to judge the astrometric accuracy, especially since 
all sources presumably are members of the star-forming region. 
The general astrometric accuracy of 2MASS is better than 071. 
A neighbouring source to YLW 15, some 8" away is fitted to 
better than 071 while the remaining 2MASS sources, all quite 
close to the limits of the field of view though in different direc- 
tions from YLW 15, are fitted to better than 072. The astrometric 
precision of the infrared data is thus estimated to be better than 
072, or at least 073. Fig. [4] shows the vicinity of the detected 
infrared source t ogether with the positions from prev ious proper 
motion analysis (ICuriel et al.l2003l ; lGirart et al.l20 04, see above) 
and the position of the weak VLBI source. The infrared emission 
appears to be unrelated to the VLBI source. 

However, the observed near infrared (1-2 /mi) emission from 
deeply embedded objects often is not directly due to the pro- 
tostar, but rather a combination of stellar light scattered at cir- 



J. Forbrich et al.: Searching for coronal radio emission from protostars using Very-Long-Baseline Interferometry 



7 



YLW 15 (VLT-ISAAC J) 



49.5 - 




Right ascension 



Fig. 4. Archival VLT-ISAAC /-band data of YLW 15 (grey-scale 
image), obtained only 1.5 months after our VLBI run (16 June 
2005). Only VLA2 is supposedly seen in this band. The VLBI 
source (marked by a circles with radii of O'.'l, U.'2, and 0'.'3; see 
text for discussion of astrometric accuracy of the infrared data) 
is apparently unrelated to the infrared emission. For reference, 
the p ositions from previous proper motion analysis are shown as 
well dCuriel et al.ll2003HGirart et al.ll2004 see text). 

cumstellar material as well a s line emission from shock heated 
regions in jets and outflows ( Wei gelt et al.l l2006). Thus, the po- 
sition in the /-band image might be just light scattered off the 
dust. 

For all sources discussed here, it follows from our observa- 
tions that the main radio emission observed with the VLA does 
not come from sub-AU scales around these protostellar sources. 
This can be explained by radio emission due to an optically thick 
thermal emission region which is resolved out at VLBI resolu- 
tion. In fact, while there are some indications of partially non- 
thermal radio emission in our target sources, there are also clear 
signs of thermal emission (see source descriptions above). 

5. Summary 

We present the results of a search for non-thermal radio emission 
on sub-AU scales towards a specially selected sample of proto- 
stars showing signs of potentially non-thermal centimetric ra- 
dio emission. The observations, partly carried out with the High 
Sensitivity Array, involving most of the largest radio telescopes 
in the world, at an aggregate data rate of 256 Mbit s _I , are the 
most sensitive yet towards this class of objects. While all sources 
in this sample are detected in VLA-only observations, we fail to 
detect compact emission in VLBI observatons except for a per- 
sistent detection of compact radio emission towards the putative 
position of YLW 15-VLA2 with a significance of 7 . . . 9<x de- 
pending on the imaging parameters. In this case, the size of the 
corona can be constrained to < 0.4 x 0.1 AU. Given the appar- 
ently fast orbital movement of YLW 15 VLA-2 around VLA-1, 
more observations will soon enable a much better determination 
of the orbit. A plausible explanation for the non-detections of 
the other sources is that thermal free-free emission from ionized 
regions (e.g., due to stellar winds or jets) is seen in the VLA 



observations which is resolved out in our VLBI data. It seems 
that the only means to overcome this problem is to look for tar- 
get sources which have favourable lines of sight allowing for a 
more direct view of the non-thermal emission in order to reduce 
the role of the optically thick free-free emission, e.g., by ensur- 
ing the clear detection of non-thermal radio emission already in 
VLA data. Since VLBI observations offer the possibility to study 
the active coronae and thereby the magnetic fields of protostars 
at high angular resolution down to < 0.1 AU, we hope that these 
observations will be possible for other candidate sources in the 
future. 
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Table 4. VLBI results 



Source 


Distance 


synth. beam (FWHM) 




map rms 1 


YLW 15 


130 pc 


2.9 mas x 0.9 mas, PA -3.6° 


0.4 AUx 0.1 AU 


15.4 ^Jy 


EC 95 


310 pc 


2.6 mas x 0.7 mas, PA -7.7° 


0.8 AU x 0.2 AU 


14.3 /Jy 


HLTau 


140 pc 


1.3 mas x 0.7 mas, PA -6.3° 


0.2 AUx 0.1 AU 


ll.OyuJy 


LDN1551 IRS 5 


140 pc 


1.3 mas x 0.7 mas, PA -6.1° 


0.2 AUx 0.1 AU 


11.4 //Jy 



1 determined from a Stokes-/ clean map with 072 side length (1024 x 0.2 mas) 
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